Two Qubits Tavis-Cummings Model Beyond the Rotating Wave Approximation: 

Degenerate Regime 
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We study the dynamics of two qubits interacting with a single mode of a harmonic oscillator 
beyond the rotating wave approximation in the ideally degenerate regime. Exact analytic expressions 
are obtained for state properties of interest, including qubit entanglement for a certain class of initial 
states of the oscillator and the qubits. Qualitative differences and similarities in the evolution of the 
qubits in the degenerate regime when the oscillator is treated quantum mechanically and classically 
are discussed. 

PACS numbers: 



I. INTRODUCTION 

Two level systems that interact with a harmonic os- 
cillator can model many physical phenomena, such as 
atoms interacting with an electromagnetic field jl|, Q], 
electrons coupled to a phonon mode of a crystal lattice 
[3], or super-conducting qubits interacting with a nano- 
mechanical resonator [4, 5], a transmission line resonator 
[1, 01 J LC circuit [8,, ^] . In the area of quantum op- 

tics, the model describing the interaction of a single qubit 
with a harmonic oscillator, called the Jaynes-Cummings 
model (JCM) ^ is extensively studied in the literature 
[To| . A generalization of the JCM to include many qubits 
not interacting with each other, but interacting with a 
single mode of a quantized harmonic oscillator, called the 
Tavis-Cummings model (TCM) [H!, [13, is also of great 
interest as it provides a framework to study collective 
properties of the qubits. 

In physical situations where the qubits are nearly res- 
onant with the oscillator and the coupling between the 
qubits and the oscillator is weak, it is a good approxima- 
tion to drop certain terms, called the counter-rotating 
terms, from the Hamiltonian describing the evolution of 
the qubits and the oscillator. Under this approximation, 
called the rotating wave approximation (RWA), the dy- 
namics of the system can be solved analytically [l[ • 

In optical setups, the near resonance and the weak cou- 
pling conditions are usually satisfied under appropriate 
conditions, justifying the use of the RWA to describe the 
dynamics of the system With the recent advance- 
ments in the area of circuit QED, it is now possible to 
engineer systems for which the near resonance and / or 
the weak coupling conditions are not satisfied jl3l4l5l |. 
To understand the dynamics under such conditions, it is 
necessary to study the evolution of the system beyond the 
RWA. Analytical approximations and numerical methods 
have been used to understand both qualitative and quan- 
titative features of the dynamics beyond the RWA for the 



tnes-Cummings [l6l - l29| and the Tavis-Cummings 11- 
model. 

In this report, we study quantum collective dynam- 
ics, in particular the entanglement properties of the two 
qubits TCM, beyond the RWA, in the regime when both 
the qubits are far from resonance with the oscillator. 
We extend recent treatments [s^ - iSTj of the role of the 
counter-rotating terms in the dynamics of entanglement 
between the qubits. 

In particular, we examine the idealized case of zero- 
gap qubits, focusing on the entanglement dynamics of 
the four Bell states with the oscillator initially in either 
a thermal state, a coherent state or a number state. A 
study of the dynamics of a degenerate qubit interacting 
with a classical field, and a discussion of a physical system 
which can be treated as a degenerate qubit, has been 
given by Shakov and McGuire |38| . 

We find the eigenvalues and eigenfunctions of the 
Hamiltonian of a two-qubit TCM in the degenerate 
regime in section |TT1 In section jilll using the Glauber- 
Sudarshan P representation [S^H^] to describe the initial 
state of the oscillator, we find the evolution of any initial 
product state of the qubits and the oscillator. The entan- 
glement properties between the qubits for a certain class 
of initial states are studied in section ITVl In section |Vl it 
is shown that the dynamics of the qubits is qualitatively 
different for the quantum mechanical and the classical 
descriptions of the oscillator and we conclude in section 



II. THE MODEL 

The Hamiltonian describing the dynamics of two 
qubits interacting with a single harmonic oscillator is 



given by [11 
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(1) 

The two qubits are assumed to be similar but fundamen- 
tally distinguishable. The transition frequency of each 
qubit (assumed to be the same for simplicity) is ujq and 
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FIG. 1. Energy level diagrams for: (a) ljo ~ ^ and (b) cjq ^ 
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the oscillator frequency is w. The coupling constant be- 
tween oscillator and both the qubits (again assumed to be 
the same for mathematical convenience) is A. The rais- 
ing and lowering operators for the harmonic oscillator are 
conventional. The qubit operators ai*'' and ai^^ are the 
usual Pauli matrices in the Hilbert space of the j*'' qubit. 
The interaction part of the Hamiltonian that couples the 
qubits to the oscillator is Hj = hX(a + 

While discussing the RWA, it is suggestive to write the 
interaction Hamiltonian as: 

Hi = hX{a + a^){a'^^^ +a^^^) + nX{a+a^){a^_^^ (2) 

If A ^ cj, cjQ and ~ one can neglect the counter ro- 
tating terms: o^ct^-* and aa^\ This is called the rotating 
wave approximation [H, The condition wq ~ w, when 
the RWA is valid is shown in Fig. [Ija). The zero-gap 
degenerate regime that we are interested in takes ljq — 0, 
an idealization of the case when uq <^ to and ^ A. 
The condition <C w is shown in Fig. [Ub). 

When LdQ — 0, the non-RWA Hamiltonian takes the 
form 



Ha = huja)a + h\{a + a)){d'^J-'^ + 4^)). 



(3) 



In order to study the dynamics of the system in the de- 
generate regime, we first find the eigenstates, |$), and 
eigenvalues, E, of Hq: 



nwa^a + h\{a + a)){a^^^ +df'^) \ |$) = E\^). (4) 



Eigenstates |<i>) will be of the form |j, m)\(t>m) where |j, m) 
are the eigenstates of (o-i^'' +d^jP) and \4>m) are the oscil- 
lator eigenstates found from Hq by replacing {cr'^x^ 
by the eigenvalue corresponding to | j, m) [24;] . 



The four eigenstates of ((ri^' -I- a)c') are: 
|j,m) = |1,±1), |1,0) and |0,0), 



(5) 



with eigenvalues 2m. In terms of the simultaneous eigen- 
states of and cf^\ cfx'^\±.) = ±|±), the states \i,m) 



can be written as: 
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(6) 



Having found \j,'m), let us now find \4>m) that satisfy 
the eigenvalue equation: 



El 



(7) 



We denote X/oj by /3 and mX/oj by /3m. Taking /3m to be 
real, and completing the square in equation JT]), we get: 



(8) 



Using the displacement operator, D{x) = exp[x(a^ — a)] 
(for real x), we can write the expression on the left hand 
of §1 as: 

(at + 2/3m){a + 2/3mMm) = D{~2/3mpaD{2/3mMm). 

As the operator D{—2f]rn)a^ a,D{2/3m) is the number op- 
erator for a displaced harmonic oscillator with eigenstates 
D{-2/3m)\N), we get from equations dH) and (O that 



\<j>m) = Dh2/3m)\N) = \N„ 



(10) 



Here |iV) is the number state of the undisplaced oscilla- 
tor, i.e. ata|iV) = N\N). 

Corresponding to the state \Nm), the energy, E, in 
equation ([7]) takes value: 



(11) 



In the discussion of two level systems interacting with a 
harmonic oscillator beyond the RWA, the use of a dis- 
placed harmonic oscillator basis was first used by Crisp 
[2(j |. The displaced oscillator states have the properties: 



{Nm\Mm') + 0. 



(12) 



The non-orthogonality condition, {Nm\^lm') 7^ 0, plays 
an important role when one discusses the spectrum of the 
system for 7^ 0. 



III. EVOLUTION OF ANY INITIAL PRODUCT 
STATE OF OSCILLATOR AND QUBITS 

Let the initial state of the system be given by the den- 
sity matrix /5(0) = Q (E) F, where Q is the joint density 
matrix of both the qubits and F is the density matrix 
of the oscillator. Writing Q in the basis \j,m), and the 
oscillator density matrix in the diagonal coherent state 
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basis (Glauber-Sudarshan P representation), p{0) takes 
the form: 

Pi*^) =^Q]m;frn'\j,m){f,m'\(g} / aP {a)\a) {al 



(13) 

where Qjm-j'm' — 0, 'tt-IQIj', m'). Evolving p(0) in time, 
we get: 

pit) = e-'"°'^''p{0)e'"°'/''\ 
= ^ J (faP{a)Qj„,,j^„,> 

xe-'^°'/^\j,m)\a){a\{j',m'\e'"°'^''\ (14) 
From Appendix A, we get 



e-'^°*/''|j,m)|a) = |j,m)|(a + 2/3^)6 



2/3„ 



where 



X exp(-2i/3™/(wt,a)) 

X exp (4i/3^(wi - sinwt)), (15) 



/(wi,a) =siny(a*e'"*/2^ae-*"*/2)^ (^g) 



Using equation (IT5|) . p(i) takes the form: 
X b',rn)(j',m'| 

® |(a + 2/3„)e-^"* - 2/3„)((a + 2/3™0e~'"* " 2/3™'| 

X exp(-2i(/3„ - /3™')/(a;t,Q!)) 

X exp {Mpl ~ Pl,){Lut - smut)) . (17) 

The dynamics of the qubits can be fohowed by tracing 
the oscillator degrees of freedom from equation ([TT)) to 
get the reduced density matrix for the qubits: 

p,{t)^TTo{m}, 

= ^ y d^aP{a)Qjm j'm' 
X \j,m){j',m'\ 

x{{a + 2/3™0e~'"* " 2/3™'|(a + 2/3™)e-'"* - 2/3„) 

X exp{-2i{l3m- l3m')f{i^t,a)) 

X exp (4i(/3,^ - Pl,)iLot - sinc^t)). (18) 

Using the relation of scalar product between two coherent 
states, {iy2Wi) = e~l''i-''2l'/2e(''2*'^i-'^2!^r)/2^ gg^. 

Pq{t) = ^ J d^aP{a)Q.j,n-.j'm' 

X |j,m)(j',m'| 

X exp (-8(/3,„ - /3„0^ sin^ ujt/2) 
X exp(-4i(/3„ - 

X exp (4i(/3^ - /S^Ol'^i - sinwt)). (19) 

The reduced density matrix for the qubits, equation (1191) . 
will be used in the next section for studying the dynamics 
of entanglement between the two non-interacting qubits. 



It must be noted that for some oscillator states, the 
function P(a) diverges, e.g. squeezed state. In such 
cases, the integral in equation (fT9|) can be ill-defined. 
To deal with such states, instead of using the P{a) rep- 
resentation, one might find it convenient to write F as: 



F= — / d^ad'^a'\a){a\F\a'){a'\ 



(20) 



and then follow the evolution of the composite system. 



IV. DYNAMICS OF ENTANGLEMENT 
BETWEEN THE QUBITS 



Using a Hamiltonian similar to ([3]), it was shown in 41 1 
that two non-interacting, initially unentangled qubits can 
get entangled due to their interaction with a common 
bath (harmonic oscillator) initially in the ground state. 
In this section, we study the dynamics of entanglement 
between the two qubits for various different initial states 
of the oscillator and the qubits. For all the initial states 
that we study, the qubits and the oscillator are chosen 
to be initially separable, and are taken to be of the form 
of equation (IT^ . Of all the possible initial states of the 
form (Uni), we choose to analyze the states for which the 
qubits are initially maximally entangled and are in one 
of the four possible Bell states: 



|vI/±) = — (|e,5)±|5,e)), 
|$±) = -i=(|e,e)±|5,ff)), 



(21) 



where ai^^\e) = |e) and (ri^^lg) — — Iff). 

From equation ([6]) , we know that for any initial state of 
the oscillator, the states \'^-) and |$_) are eigenstates of 
Hq. So, the reduced density matrices for the qubits will 
not change if they start in the state l^*-) or |<&_) and the 
qubits will remain maximally entangled no matter what 
the initial state of the oscillator is. 

The local unitary operator, ^i*^^'^, commutes 

with the Hamiltonian and connects the states 1^*+) and 
l$+): 



0. 



(22) 



Since any measure of entanglement calculated for two 
different states that are connected by a local unitary 
transformation has the same value [42] , the entanglement 
dynamics between the qubits for the initial state |^'+) 
will be the same as entanglement dynamics between the 
qubits for the initial state |$+). One should note that 
this result is valid for any initial state of the field. Be- 
cause of this, we explicitly study only the states for which 
the initial state of the qubits is 1^*+). 
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When the qubits start in the state_[^'+), it turns out 
that the expression for concurrence (ijl, a measure that 
quantifies entanglement between two qubits, has an an- 
alytic general form that is valid for all the initial states 
of the oscillator. Wc first derive that general expression 
and then study explicitly the cases for which the initial 
state of the oscillator is either a thermal state, a coherent 
state or a number state. 

For the state |^'+), the initial density matrix for the 
qubits, Q, takes the form: 



1 

1 ( 

2(0 00 

-10 1 



IN 



(23) 



where the rows and columns are defined in the basis: 
|1,1), |1,0), |0,0) and |1,-1). We denote the initial 
state of the oscillator by F, which is characterized by 
the Glauber-Sudarshan function P. As time evolves, by 
using equation (jl9p , we calculate the reduced density ma- 
trix for the qubits in the \ j, m) basis and get: 



1 








-m 


























i{tr 











(24) 



where 



/(i) = J (faPia)exp{-8iPf{ujt,a)) 
X exp (-32^2 sin^ ('^i/2)) ■ 



(25) 



In order to calculate the concurrence in the usual form 
between the two qubits, we need to write Pq(t) in the 
|e,e), \e,g), \g,e) and \g,g) basis: 



n - u{t) iv{t) iv{t) 

-iv{t) l + u{t) l + u{t) 

-iv{t) l + u\t) l + u{t) 

\1 - u{t) iv{t) iv{t) 




(26) 



The terms u{t) and v{i) are the real and imaginary parts 
of I{t) respectively. For this density matrix, the concur- 
rence is given by: 

c{t)^\m\, 

d^aP{a) exp {-8i(3f{ujt, a)) 

X exp (-32/3^ sin^ {ut/2)), 
= P(8/3sin(a;i),8/3(cos(a;t) - 1)) 

X exp (-32/32 sin2(tjt/2)). (27) 

The function P, defined in equation (|27p . is the Fourier 
transform of P: 

fca) = / dxdyPix, y) exp {-i{kix + fcaj/)), (28) 



with X and y being respectively the real and imaginary 
parts of a in the integral (|25|) . From equation (f27l) . we 
see that knowing P(a) or its Fourier transform, P, that 
characterizes the initial state of the oscillator, we can 
follow the evolution of entanglement between the qubits. 

Qubits initially in state |^+) and oscillator in 
the thermal state: The dynamics of a system interact- 
ing with a thermal reservoir is generally associated with 
decoherence of the system. However, following the work 
of Bose et al. [44], it was shown in [i^ that entangle- 
ment (which is a manifestation of coherent superposition 
of states of a multi-partite system), can be created be- 
tween the two qubits of the Tavis-Cummings model with 
the oscillator initially prepared in the thermal state. En- 
vironment induced entanglement is an area of active re- 
search jillieli. 

Entanglement dynamics of the two-qubits Tavis- 
Cummings model with the oscillator initially in a ther- 
mal state is investigated in [47l - l50j . In particular, Roa et 
al- [53, studied the entanglement dynamics between the 
qubits in the dispersive coupling regime within the rotat- 
ing wave approximation, where the qubits were taken to 
be far off resonance from the oscillator and the coupling 
between the qubits and the oscillator taken to be weak. 

To study the entanglement dynamics between the 
qubits in the degenerate regime, beyond the RWA and 
for arbitrary coupling strengths, for the initial state of 
the qubits being and that of the oscillator being 

a thermal state using equation (127)) . we note that for a 
thermal state P{a) is given by |5l|: 



Pth{a) = 



1 



■exp(-|a|V(n))- 



(29) 



In the above equation, (n) is the average number of ex- 
citations for a given temperature which characterizes the 
thermal state. We note that as P{a) is a gaussian of 
width \/(n), its Fourier transform, P, is also a gaussian 
of width l/y/{n). 

Using Pth{a) in equation (P7)) . we get: 

Cthit) = exp (-32/32 gjj^2 (^t/2)(l + 2 (n))) . (30) 

In Fig. [21 Cthit) is plotted for various values of (n) and 
/3 = 0.05 and /3 = 0.1. We notice that entanglement 
between the qubits oscillates periodically with time pe- 
riod l-Kjuj. Within each time period, the entanglement 
between the qubits decreases more rapidly from its max- 
imum value as one increases (n) or /3. The minimum en- 
tanglement between the qubits, exp (—32/3^(1 -|- 2 (n))), 
reduces as (n) or /3 is increased. These are consequences 
of P becoming narrower as (n) or /3 increases. 

Qubits initially in state |^'+) and oscillator in a 
coherent state: The dynamics of a single qubit interact- 
ing with a coherent state of the oscillator has been stud- 
ied extensively in the literature and has revealed many 
important phenomena that are unique features of quan- 
tum mechanics such as collapse and revival of atomic 
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FIG. 2. Time dependence of entanglement between the two 
qubits in the degenerate regime given by Cth(f) for various 
value of in). Top: /3 = 0.05 and bottom: (3 = 0.1. 



coherence [521 ^ entanglement and disentanglement of the 
qubit with the oscillator and in process preparation of 
the Schrodinger Cat state etc. Within the RWA, 
the collapse and revival of a two-qubit TCM was studied 
in (5^ . and the entanglement dynamics for various dif- 
ferent initial states of the oscillator, including coherent 
state, was done in (ssj . 

For the oscillator initially in the coherent state |ao), 
P{a) is given by 51): 



Pcoh{a) = (5^(a - ao). 



(31) 



With Pcoh{c^), I{t) in equation ([25]) is evaluated to be: 



Icoh{t) = e^p{-8il3f{ujt,ao)) 
X exp 



(32) 



From equations (|27|) and (|32|) . and the fact that f{u!t, ag) 
is real, the concurrence for the coherent state case is 
found to be: 



Ccohit) = exp (-32/32 ^-^2 (^i/2)). 



(33) 



We see from equation ([33| that the entanglement between 
the qubits does not depend on the oscillator strength, ao, 
and only depends upon /3 and lu. The fact that Ccoh(t) 
does not depend upon ag is a consequence of P being 
just a phase as Pcoh{ct) is a delta function. In Fig. [3l we 
plot Ccoh{t) for various values of /3. We see from Fig. [3] 
that Ccoh{t) is periodic with time period 27r/cj. As in the 
thermal state case, the entanglement between the qubits 
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FIG. 3. Time dependence of entanglement between the two 
qubits in the degenerate regime given by Cco{i) for various 
value of /3. 



decreases more rapidly from its maximum value as one 
increases /3. 

Qubits initially in state |^'+) and oscillator in 
a number state: To study the entanglement dynamics 
between the qubits in the degenerate regime when the 
oscillator starts in the Fock state, we note that P{a) for 
a Fock state |iV) is given by [Slj: 



PN{a) 



exp [aa*) ( d 



2N 



5\a) 



(34) 



With Pn{ol) in equation ((25|) . I{t) becomes (see Ap- 
pendix B): 



lN{t) 



N 



1 



N 

^ \ml (N - m)\ 

m=0 ^ / ^ ' 

X exp (-32/3^ sm^{ujt/2)), 

((8/3 sin {Ljt/2)f) 
X exp(-32;32 sm^{ujt/2)), 



(8i/3sin {Ljt/2)f^-"'^ 



(35) 



where Ln{x) is a Laguerre polynomial. Using equations 
([35]) and ([27]) . we calculate the concurrence for the num- 
ber state case: 



CN{t) = |LAr((8/3sin(L^t/2))2)| 

X exp (-32/32sin2(wt/2)). 



(36) 



The concurrence, CN{t), is plotted in Fig. |4] for var- 
ious values of N and /3 — 0.05 and /3 = 0.1. We see 
from the figure that C]\f(t) is periodic with time period 
2tt/ui. Within each time period, for certain values of N 
and f3, the qubits momentarily become disentangled be- 
fore getting entangled again. As N 01 (3 is increases, 
entanglement between the qubits gets sharply peaked at 
times that are integer multiples of 27r/aj. 

Dependence of entanglement on coupling 
strength: As discussed in section [III for the RWA to 
be valid, one has to ensure that A ^ a;,wo. So, given uj 
and Wq, there is an upper limit of A for which RWA is 
valid. In contrast to the RWA case, the results in this re- 
port are correct for all values of A. We see from equations 
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FIG. 4. Time dependence of entanglement between the two 
qubits in the degenerate regime given by Cjv(t) for various 
value of N. Top: /3 = 0.05 and bottom: /3 = 0.1. 



FIG. 5. Entanglement between the qubits as a function of /3 
for (a) thermal state, (b) coherent state, and (c) number state 
evaluated at ujt/2 = 7r/2. 



(|33p and that the entanglement depends on the 
coupling strength and the oscillator frequency through 
the parameter sin^(a;t/2). This parameter is periodic 
with time period T — 2ii/uj. At times that are integer 
multiples of T, the qubits are maximally entangled for all 
coupling strengths. This can be seen from Figs. [2l[3]and 
m To understand the role of the coupling strength, A, in 
the entanglement between the qubits, we plot in Fig. [S] 
the concurrence for the thermal state, the coherent state 
and the number state as a function of the dimension-less 
coupling strength /? evaluated at half the time period, i.e. 
t = it/u). We have: 

Cth{^/uj) = exp (-32/32(1 + 2 (n))), (37) 
Ccoh (V^) -cxp(-32/32), (38) 
CAr(7r/w) ^\Ln I exp (-32/3^) . (39) 

We see from the above set of equations that for all 
the three initial states of the oscillator, concurrence de- 
creases exponentially as P increases. Equation (37), for 
various values of (n), is plotted in Fig. [Sj^a) from which 
we see that the concurrence decreases faster as the mean 
thermal excitation number (n) is increased. For the co- 
herent state, equation (38) is plotted in Fig. [Sfb) from 
which the exponential decrease of concurrence is evident. 
In Fig. Eljc), Cn{-k/uj) is plotted for various values of N 
where we see that the exponential decrease of concur- 
rence is modulated by the Laguerre polynomial. 



V. CLASSICAL-QUANTUM COMPARISON 

It is well established within the rotating wave ap- 
proximation that the predictions about the dynamical 
properties of the qubits interacting with a classical har- 
monic oscillator are qualitatively different from the pre- 
dictions when the oscillator is treated quantum mechan- 
ically. Contrary to the classical case involving a pre- 
scribed driving field, the quantum mechanical treatment 
for the corresponding coherent field predicts the qubit 
population inversion to show collapse and revival (s^ . 
These arise from the infinitely many, but discrete, Rabi 
frequencies for qubit evolution associated with the infi- 
nite set of Fock states defining a coherent state. The 
quantum mechanical predictions have been confirmed ex- 
perimentally [1^. In view of this, it is a question whether 
there are qualitative differences between the predictions 
that can be made on the basis of the non-RWA calcu- 
lations in Sec. IV and a set of corresponding non-RWA 
predictions to be made when the oscillator is treated as 
classical, while the qubits remain quantum mechanical. 

A mixed quantum-classical theory of this kind is usu- 
ally called semi-classical and the oscillator is taken as an 
external influence. Here we want something a bit differ- 
ent because in Sec. IV the oscillator has not been treated 
as an external influence, but rather as open to back reac- 
tion from the qubits. Jaynes flrst addressed this issue of 
back reaction in the cavity QED context with a so-called 
"neo-classical" theory P, Is^l in which the oscillator (cav- 
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ity electric field in his case) was driven by the qubit via 
its quantum average values. In this way, although in- 
fluenced by the quantized qubits, the oscillator remains 
unquantized for all time and takes only numerical, not 
operator, values. 

We will adopt the Jaynes prescription in order to in- 
clude back reaction. The simplest approach is to derive 
equations of motion for the various dynamical quantities 
using the original Hamiltonian ([3]). For the oscillator, in- 
stead of working with non-hermitian operators a and , 
we work with the real dynamical quantities, momentum 
p and co-ordinate q defined in Appendix C. We then re- 
place operators by numerical variables, which we indicate 
simply by removing the carets. These variables evolve 
from classically specifiable initial values p(0) and (?(0), as 
driven by expectation values of the qubit cTx variables: 



q=p + h\^2inw{a^^' + a. 

P = 



(2)\ 



0, 



(40) 
(41) 
(42) 
(43) 
(44) 



As it happens, the states we have worked with have 
zero expectation of the operators, which don't evolve 
themselves, so the oscillator evolves completely inde- 
pendently. For the qubits starting in the state 
there are only four non-zero density matrix elements: 
(l,±l|p,|l,±l) and (l,±l|/5g|l,Tl)- The diagonal den- 
sity matrix elements do not change in time because 



CTa;^^ + o'x^' = 0. To find the evolution of the off-diagonal 
elements, we note that: 

(l,l|p,(t)|l,-l) =Tr{p,(t)|l,-l)(l,l|} 



^Tr{p,{0)Si'\t)S':>{t)}, 



where 



(45) 



(46) 



From equations (|43)) and (|44| , we get the following equa- 
tion for the required expectation values 

l{si'\t)S^_^\t))=-4^X^p{Si'Ht)S^_^Ht)). (47) 
Solving equations (gO]), (HI]) and (gT]) we get 

X exp sin {ujt/2) ] 

X exp (^-4il3y/2/hLu p{0) sin(a;t)) ,(48) 

where g(0) and p{0) are the initial position and momen- 
tum of the oscillator respectively. For the initial state 



of the qubits being \<f+), {S':^\o)S^^\o)) = -1/2. The 
entanglement between the qubits is then given by (equa- 
tions (HZl) and (gS])): 



C{t) 



2\{si'\t)S^^\t))\. 



(49) 



Since the time dependent terms in equation ()48p enter 
only as phases, the entanglement between the qubits do 
not change and the qubits remain maximally entangled. 

A crucial point used in deriving equation (j48p is the 
fact that the oscillator is assumed to have well defined 
initial position and momentum. Because of the uncer- 
tainty principle, within quantum mechanics, no state can 
be a simultaneous eigenstate of position and momentum. 
In view of this, it is a question whether with a classical 
description of the oscillator one can or cannot get modu- 
lations in entanglement of the qubits if the initial position 
and momentum of the oscillator are not well defined. We 
now show that by averaging equation (1481) over appropri- 
ate joint probability distributions of ^(O) and p(0), say 
P(go,Po)i we get expressions for concurrence that match 
with the quantum mechanically derived results for ther- 
mal states (equation ([50)) ) and the coherent states (equa- 
tion (1331)). 

The classical probability distribution P(go,Po) of a 
harmonic oscillator in thermal equilibrium is gaussian. 
For an appropriate probability distribution of the oscil- 
lator in a coherent state, lag), we use the Wigner distri- 
bution of the coherent state pV\ which is also a gaussian. 
So, for the purpose of thermal states and coherent states, 
we take P{qo,po) to be of the general form: 



1 



nAqAp 



exp 



(go ~ go)^ _ {pq -Pof 
(Ag)2 {Apf 



m 

where Ag and Ap are the initial uncertainties in position 
and momentum, centered at go and po respectively. Tak- 
ing the absolute value after averaging equation (g5)) over 
P{qo,Po), we get 



C{t) = 2 



dqodpoP{qo,Po){S^+\t)S^_^\t)) 



io{Aqy^ 



= exp 1^-32/32 sin'' {iot/2) 
X exp ( -8/3^ sin^ (ujt) 



h 

{Apf 



(51) 



From the above equation, we see that the concurrence 
does not depend on the average values go and po and 
only depends on the uncertainties, Ag and Ap. If both 
the uncertainties, Ag and Ap, are not simultaneously 
zero, the concurrence does change in time. 

For a thermal state characterized by temperature T, 
Aqth = ^/2KTJ^ and Apth = V2KT where K is the 
Boltzmann constant. Using this in equation (I5ip . we get: 



KT\ 

Cth (t) = exp ( -64/32 gij^2 (^^/2) — j . (52) 
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In order to compare equation (j52p with the quantum re- 
sult (eauation ([5n)) ). we note that in the high temperature 
hmit, i.e. when KT 3> fiw, the average excitation num- 
ber (n) = KT/huj. Now, replacing the factor (1 -I- 2{n)) 
in the exponential of equation (pO)) by 2KT/fiu!, we get 
exactly the same expression for concurrence as derived 
with the classical description of the oscillator, equation 

o. 

A coherent state is a minimum uncertainty state for 
which Aqcoh = \/W^ ^Pcoh — \/^- Putting this 
in equation ([5T|). we get: 

Ccohit) = exp (-32/?2 sin^ {ujt/2)), (53) 

which exactly matches with the corresponding quantum 
result, equation p3p . 

The above results show that the modulations in the en- 
tanglement between the qubits is a consequence of initial 
uncertainty in qo and po, and not an intrisically quantum 
property. However, the uncertainty in the simultaneous 
assignment of values to canonically conjugate variables is 
inherent in quantum mechanics and cannot be avoided. 
This fundamental uncertainty is reflected in the general 
equation for concurrence that is valid for all quantum 
states, equation (l27l) . from which we see that for all states 
there is a time dependent factor, exp (—32/3^ sin^ (ajt/2)) , 
which corresponds to the modulation of entanglement for 
the minimum uncertainty state, i.e. the coherent state. 

VI. CONCLUSION 

In this report, we studied the dynamics of the two- 
qubit Tavis-Cummings model in the degenerate regime 
beyond the RWA, obtaining exact analytic expressions 
for state properties of interest, including qubit entangle- 
ment. The general form of the time evolved density ma- 
trix of the composite system for any initial product state 
of the qubits and the oscillator was derived. The qubits 
initially prepared in the |^-) and 1$-) states were shown 
to remain maximally entangled throughout the dynam- 
ical evolution. The dynamics of entanglement between 
the qubits initially prepared in the |vE'+) state was shown 
to be the same as that for the qubits initially prepared 
in the state |$+) for any initial state of the oscillator. A 
general expression for concurrence, quantifying the en- 
tanglement between the qubits initially prepared in the 
state |^+) and the oscillator in any state characterized 
by the Glauber-Sudarshan function P{a) was derived. 

When the initial state of the oscillator is either the 
thermal state, coherent state or number state, the entan- 
glement dynamics between the qubits was shown to be 
periodic with time period 2t: /lu. This periodic behavior 
is a consequence of the fact that our system consists of 
displaced harmonic oscillators with frequency oj. A gen- 
eral feature of entanglement for all initial states of the 
oscillator was a more sharply peaked entanglement as the 
dimension-less coupling strength /3 increased. This fea- 
ture can be understood as a manifestation of the common 



factor exp (— 32/3^sin^ {ujt/2)) in equation (|27l) . which is 
present for all initial states of the oscillator. 

When the initial state of the oscillator is a coherent 
state, it was shown that the entanglement between the 
qubits does not depend on the oscillator strength uq . The 
periodicity of entanglement between the qubits with a 
period independent of the strength of the coherent state 
was reported in where the dynamics between the 
qubits and the oscillator was studied in the dispersive 
regime within the RWA. Although in \5d\ the periodicity 
was found not to depend on the strength of the coherent 
state, the amplitude of oscillations of entanglement de- 
pended on ao in contrast to the results found here in the 
degenerate regime. 

In section |Vl we studied the dynamics of the state of 
the qubits by treating the oscillator classically. The back 
action of the qubits on the oscillator evolution equations 
was included in the spirit of Jaynes' original neoclassical 
treatment |57j] . It was shown that temporal entanglement 
modulation is also present classically, but only if there is 
uncertainty in the initial position and/or momentum of 
the oscillator. 

In all the calculations in the current paper, the qubits 
were assumed to be perfectly degenerate. This is obvi- 
ously an idealization. A more realistic scenario would 
have quasi-degenerate qubits with cjo 7^ but still sat- 
isfying the condition ujq <^ oj, A. The existence of com- 
parable analytic expressions for the collective properties 
of such qubits with non-zero ujq is an important question 
for future work. 
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APPENDIX A: EVOLUTION TERM 

= |j,m)i3(-2/3™)e-^'^*^'"73(2/3„)|a)e4''5""*, 
= |j,TO)|(a + 2/3„0e-*'^*-2/3„) 

X exp(-2i/3„siny(a*e^'^*/2 -^ae-*"*/^)) 

X exp (4i/3^(iLjt - sinujt)). (54) 



APPENDIX B: CALCULATION OF lN{t) 
We define 

ij^{t) = /at (i) exp (-32/32 sin^ (a;i/2)), (55) 
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and evaluate lN{t) for -PAr(a) given by equation (|34|) : 
lN{t) = J d^aPNia) 

X exp(-8i/3siny(a*e'"*/2 _^ae-*"*/2)) 

If ( d"^^ 

= _y d^«exp(|an(^^-,^^<S^(a) 

X exp(-8i/3siny(a*e*'^*/2 ^ae-*"*/2)) 



d^aexpdap) ^^^<s2(a; 



X exp {A{a*a + a6)). (56) 



where 



2 



a = e*"*/2^ and 6 = e-'"*/^. 
Using the identity 



(57) 



d^a 



2N 



d 



2N 



/(a, a*) |q=q.=o, 



(58) 



we get 



1 a- 



,2N 



]a\\A{a* a+ab) 



AM^aW^ "'^ 'J 1—0, 



1 9 



TV 



[(a*+A6)^e^""*]U.^o 



AT! da*^ 

N 



Nl ^-^ \m J \ da 

m=0 \ / \ 



/ QN — m 



V da' 



N 



m=0 \ / ^ ^ 

i-fn, r-8./3sin^)^(^--). 
TV! \mj {N-m)r ^ 2' 



(59) 



APPENDIX C: HAMILTONIAN CONTEXT 

For an exphcit example of context for the generic 
Hamiltonian (|3]) we provide a version of the original treat- 
ment by Jaynes (57| . In this case we are describing a ver- 
sion of cavity or circuit QED. The electromagnetic field 
mode is interpreted as an oscillator, and the associated 
field Hamiltonian is written in terms of real dynamical 
variables, momentum p and co-ordinate q as: 



H 



-E 

2 ^ 



2 fl2 ) 



(60) 



where and qm are the expansion coefficients of the 
vector potential and the transverse electric field: 



A{r,t) = - \/eo E (^m )gm [t) , 

)Pm(i), (61) 



or their equivalents in a resonant circuit. In QED these 
are operators obeying canonical commutation relations. 
The function Um(fm) is the orthonormal vector func- 
tion (i.e., incorporating polarization) that defines the m*'' 
mode of the field, which has frequency a;„i. 

The interaction Hamiltonian is derived in the familiar 
way from the dipole -field coupling: Hint — —d.E. For a 
single excited mode labeled a, and a qubit with transition 
dipole matrix element di2, this takes the explicit form 



Hint = -{dl2-Ua{0))^/£0Pa{t)^x, 

= hX{d + d'')ax, 



(62) 



where we have introduced the coupling constant A used 
throughout, as well as the oscillator amplitudes via Pa — 
y^haj/2{d + at) and ujaqa = i\/^al'2(a — a^). As usual, 
the Pauli matrix expresses the operator character of the 
two- level dipole moment. In Sec|Vl we removed the hats 
and regarded Pa, (Jq, a, and as classical numerical vari- 
ables. 
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